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^Nj ^ In this paper, we study exactly the thermal conductance for a low dimensional system repre- 

sented by two coupled massive Brownian particles, both directly and via a Green-Kubo expression. 
Both approaches give exactly the same result. We also obtain exactly the steady state probability 
distribution for that system by means of time-averaging. 
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(U ; I. INTRODUCTION 

Exact results, in physics, play an important and useful role as a reference for other methods. For instance, they 

■ can be used to study specific features of models which are not easily accessible to approximative methods, such as 
^ \ computer simulations. However, non-trivial exact results are few and difficult to come by: in the literature of transport 

: phenomena there are not many exact calculations for transport coefficients based on the mechanical parameters of 
the systems under observation P, 0, d, 0, H B S B @- Indeed, the calculation of transport coefficients is one of the 
most important goals of non-equilibrium physics. 
I In fact, the rigorous derivation of Fourier's Law for bulk Hamiltonian systems is still under debate The soluble 

harmonic models used to evaluate the thermal conductance do not reproduce the necessary scattering of the energy 
unless local thermal reservoirs (which act in part as effective scatterers for the energy incoming on them) are coupled 
Q \ to the bulk sites. Other models, such as mass disordered ones have been proposed but there is evidence that mass 
disorder alone will not give rise to Fourier's Law in 2D ^1!]. The presence of anharmonicity on the coupling level would 
probably be sufficient for energy to scatter and diffuse from site to neighboring site. However, the technical details 
for obtaining exact rigorous results are, as far as we know, too difficult to overcome at the present level 0]- 
(y-^ On the other hand, for Equilibrium Statistical Mechanics the probability distribution for a given system can be 

00 found, given that it obeys Liouville's theorem [13] and the external macroscopic constraints [ll|, by means of an 
' ensemble of points in phase space, when that system is ergodic. However, ergodicity is not a necessary condition for 
I obtaining the long times the stationary probability distribution since it is always possible, at least in principle, to 
\^ , obtain the time-average for any physical quantity during the realization of an actual experiment. 

What distinguishes time averaging from other exact methods derived from the solution of a Fokker-Planck formalism 
On . is that time averaging can take into account, exactly, all the orders of the moments of the dynamical variables. The 

■ Fokker-Planck formalism is exact only to the second order moments while time average is akin to the stationary 
^ solution of the Kramers-Moyal equation correct at all orders of moments [l^ . For simplicity sake we use Gaussian 

r white noise in the present work but the method is readily generalizable to any type of noise, given that all its moments 
r> \ are known. 

■ Our present goal is to study exactly the validity, and consistency, of some methods used in the derivation of 
" " ' transport coefficients (namely the thermal conductance) for small classical systems. The reasons for using small 

systems in our model are manifold. Firstly, exact calculations become feasible. Secondly, macroscopic transport 
properties associated with large systems must have a mechanical counterpart in small systems. Macroscopic flows of 
mass, momentum or energy are the effect of the averaging of the action of microscopic forces, and work. Thirdly, small 
systems are interesting per se. There are difficulties inherent to small systems due to the fact that one cannot take 
the thermodynamic limit that averages out many problems associated with solving the dynamics of large systems. 
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similarly to the Law of Large Numbers that arises from the summation of many random variables ■ We can also 
take into account exactly the effects of the inertia and of a possible non-Markovian nature for the noise. 

For deriving the exact time-averaged thermal conductance, we choose a simple system which is capable of non- 
trivially transmitting heat between its constitutive parts, via mechanical work, when subjected to a gradient of (white 
noise) temperature: a system of two coupled Brownian particles (BP) with the coupling constant k acting as the sole 
information channel between the particles. In the present case, the system dynamics is linear and the linear response 
treatment shall be proven to be exact, as we will see in the following. 

We can apply a method previously used by the authors [3, to the study of thermal conduction between 
the coupled BP system. The exact time-averaging method of References [13, [S] is capable of obtaining exactly 
the stationary probability distribution for Brownian particles submitted to white and colored noise. In particular, 
by submitting a single massive Brownian particle (BP) to two different thermal contacts, at distinct temperatures 
(similarly to some glasses that are subjected to thermal vibrations and structural modifications represented by distinct 
noise functions at different time-scales (l6j). we can keep it from reaching thermal equilibrium [isj . 

Since the present model is effectively zero dimensional, the thermal conductance between the Brownian particles is 
defined simply as the energy flow per unit time per temperature difference between the particles, i.e., the conductive 
flow of energy (for particle 1, from particle 2) is defined as ji^2 = {Ti^2 — ?2,i), where k is the inter-particle thermal 
conductance in first-order approximation. Indeed, there has been some recent developments in treating finite systems 
that can be adapted to the problem under study [TtI . [18| . In that case, the transport coefficient k is obtained via a 
convenient Green-Kubo formulation. The calculation of the transport coefficients by this method can be an interesting 
starting point for the study of more complex models, such as polymers subjected to gradients of temperature [li- 
lt will also provide an important test for the choice of flow variable appropriate for such models. However, in order 
to avoid the rather artificial construction of an ensemble of reservoirs that need to be coupled to the particles along 
the linear polymer (harmonic crystal), a generalization of the method will be needed to include non-linearities on the 
potential. This way, a much more realistic picture of thermal conduction will be obtained. 

The thermal conductance between two particles is not a well defined macroscopic quantity since we are far from 
the thermodynamic limit and cannot define a macroscopic (and diffusive) flow of heat. However, it is clear that if a 
(classic) macroscopic system is partitioned into two parts, energy conduction is realized by the interactions (work) at 
the interface. 

Furthermore, we add a periodic variation of the temperature of the Brownian particles. This is an interesting 
effect that can lead to the appearance of currents for systems presenting asymmetries in the potential energy 2^ ■ 
Periodic oscillations of different types are capable of creating currents [2lj, [22| in the case of zero average forces acting 
on the particles. The combination of ratchet-type potential energy and periodic time oscillation for the temperature 
has been extensively studied [13, [l^l . 

This paper is organized as follows: In Section II we define the model. In Section III we explain the method of 
time-averaging and show the main contributions to the probability distribution. In Section IV we calculate the time- 
averaged steady-state distribution for the non-equilibrium conditions. In Section V we obtain the thermal conductance 
and in Section VI we discuss our main conclusions. 



II. EXACTLY SOLVABLE MODEL 



Our model consists of two massive Brownian particles (BP) coupled by an harmonic potential and subjected to 
white noise at distinct temperatures. This could be interpreted as two atoms in a crystal, coupled by a harmonic 
potential. 

Despite the reduced number of variables, the present system contains the main ingredients of more complex models. 
In it, we can define the energy transfer as the microscopical work, that in macroscopical systems become the internally 
transferred heat. In the following, we describe the model in detail and, using time- average techniques [ij, [l^, we 
calculate exactly the probability distribution for the relevant Brownian variables. 



A. Langevin-type equation 

The system composed by two coupled punctual and massive BPs is described by the equations: 

Xa{t) = Va{t), (1) 
maVa{t) = -k{Xa{t) - X,3{t)) - k' Xa{t) --faVa{t) +ria{t). (2) 

Gaussian behavior is to be expected for the probability distribution for the time-averaged stationary state, according 
to previous published works (in special see sections 1.3.E.2 and 2.2.E.2 in Ref. [24|). 



3 



For simplicity, we make: toi = m2 = m, k" = k' , and 71 = 72 = 7- Thus the equations can be written as: 

mXa(t) = -k (Xa{t) - Xpit)) - k' Xa{t) - JXait) + TJait), (3) 

where a, /3 = 1, 2, a /3, and the initial conditions are: 

a;i(0) = a;2(0) = t;i(0)=«2(0)=0. 



B. Noise properties 

Both white Gaussian noise terms can be defined in terms of their two lowest two cumulants: 

M)) = 0, 

where the modulated temperatures above are given by 

Ta{t)^fo, + sinKO]', 



(4) 
(5) 

(6) 



for a = 1,2 and | | < 1 ■ 

The oscillating temperatures, in other models, can induce very interesting effects such as sendingheat fluxes against 
gradients of temperature [1^ or directed fluxes of particles in periodic potentials [l^ [lO, [1^ [IB, US] ■ 



C. Laplace transformations 

Taking the Laplace transformations of Eqs.JT]) and ^ yields 



(ms^ + 7s + fc + fc') xi(s) = k X2{s) + f}i{s) , 

vi{s) = sxi{s), 
(ms^ + 7S + fc + fc') £2(5) = kxi{s) + fi2{s), 

V2is) = SX2{s). 

Defining T{s) = ms'^+'js + k + k' and rearranging Eqs.(l7]) to (fTO|) . one finds that: 



ii(s) 
vi{s) 
i2(s) 

W2(s) 



A(s) ?7i(s) + A(s)ry2(s), 
sA(s)^i(s) + sA(s)fy2(s), 
A{s)m{s)+A{s) 7)2(5), 
sA(s) ?7i(s) + sA(s)fy2(s), 



where: 



A(s) ^ 
A(s) ^ 



r(s) 



r2(s) -fc2' 

k 



r2(s) -fc2' 

The Laplace transformation for the independent noise variables is given by (a = 1,2): 



1 



2 A^ujct 



i{qi + gj) + 2e + g^) + 2e]2 + ^2 



2v4^c.2 



+ qj) + 2e] + g,) + 2e]2 + 4c^2; 

All integration paths are the same and shown in Fig[TJ 



(7) 
(8) 
(9) 
(10) 



(11) 
(12) 
(13) 
(14) 



(15) 
(16) 



(17) 
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FIG. 1: Integration path for the equilibrium distribution, Eq. (|T9 



III. PROBABILITY DISTRIBUTION 



Here we show some of tlie steps (more detail can be found in references [ij, [l^) to obtain the expression for the 
instantaneous probability distribution for the system of coupled Brownian particles. 

Time averaging, for a supposedly convergent distribution, is defined and calculated as in Refs. [l^. [Tsj: 



/ ^ J™ o / dtf{t)= lim 2 



fit) 



From the definition for the instantaneous probability distribution: 

p{xi,Vi,X2,V2,t) = {S{xi ~Xi{t))S{vi ~Vi{t))S{x2 -X2{t))5{v2 - V2{t))) , 

it is possible to show that [11] : 



p''''{xi,Vi,X2,V2) 



= lim lim 

oo 

X E 

/,m,n,o— 

f + OO ' 

X 



dQi dQ2 dPi C?^P2 ^tQlal+^Pl^l+^Q23:2+^P2f2 

, 27r 27r 27r 27r 
i-iQiY (-zPi)™ {-iQ2r {-^P2Y 



n 



n dqif -TT dpih 
1-K 27r 



h=\ 



-Q dq2i -Q Ap2k_ ^ 



z 



^ 27r J--'- 27r 

°° j=i fe=i 

rn+o+l 



ELi(*9ia + + Ebli(«Pifc + e) + Ec=i(*92c + e) + ELi(*P2<i + e) 
/=i h=l j=i 



k=l 



(18) 



(19) 



where the integration paths for the {q,p) variables are given in Fig[TJ 

An interesting case is the study of the average {xi{t)). The time periodicity of the noise [Ti^2(t + T) = Ti^2{t)] gets 
translated into a periodicity of distribution [p{x,t + T) = p(x,t)], and by consequence also of the averages for the 
variables, e.g. (^(t + T)) = (^{t))- This is consistent with a Fokker-Planck treatment of the distribution found in 
the literature 3 13, H^, [H, [13] . We now show that the instantaneous distribution described by Eq. is indeed 
consistent with periodicity in time. 

From reference [l^], it is clear that {x^{t)) is periodic in time with period T: 

{x\{t + T))— j dxip{xi,t + T) xl — J dxip{xi,t) xl — {xl{t)}, 
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where the periodicity of the probabihty distribution directly imphes that of the average. We assume that T2 = 0, 
so the present model and the one from ^] coincide. 

The average (for long time, after the memory of the initial conditions has already faded out) reads: 

{x\{t)) = / dxip{xi,t)x\ 



$i^e*[^«^+*«^+2^I(ii(*gi + e)i,{iq2 + e)) 



— lim , 

e^o7_^ 27r 27r 

By integrating the last equation above over the poles of {f]ifi2) we observe that apart form the term of order 0(^4"), 
the terms carrying the contribution from the sine bring a dependence such that iqi + iq2 + 2e — Huj. By obtaining 
the residue, we are left with exponential terms of the form e^'"*, which are periodic in time with period T. 

An interesting aspect of periodically varying noise is that it modulates the time behavior of distributions and 
averages, i.e., there are no stationary constant values for any of the moments of the Brownian variables: the moments 
are periodic functions of time. In this case, the time average we use corresponds to averages of these quantities taken 
over a period of the noise for very long observation times. In general lines, the next steps consist into expressing the 
Laplace transforms for the dynamical variables (x, v) into functions of the averages of the Laplace transforms of the 
noise. 



IV. TIME-AVERAGED STEADY STATE DISTRIBUTION 



A. Contributing terms 

The main contributions for the probability shown in Eq. (1191) comes from a typical integration is of the form p^ . \li 



2 TT 2 TT z — i{qi + qj + 2e + I 



where 



To understand what causes a term in Eq. (|19p to contribute to the time-averaged steady state probability distribution 
it is necessary to observe that for a typical integration, such as the above term, there is a factor I{z) 

I{z) 



z - i{qi + q,j + 2e + 0) ' 

that, in the limit z — > 0, will vanish if there is the presence of any finite terms on its denominator (due to the residue 
calculations around the poles of the rest of the integrand). Only integrations that eliminate all the pairs of q's in the 
denominator of I{z) will transform it into I{z) = z/ z = 1. This is a necessary condition for any of the integrations 
done below to contribute to the probability distribution. Each integration brings the corresponding multiplicative 
factors that need to be dealt with. 



We define: 



'^'^ " 47(fc'2 + 2fcA:') 4 [m fc2 + (fc + A:') -^2] ' ( 1) 

^ fc^ [m(fc + fcO+72] 

'^^^ A-fk' {2k + k')[mk^ + {k + k')-i^y ^ ' 

k 

Qaa = Qaa = , . ^o7. , (23) 



1 \ mP + 2{k + k')^^ \ 
'^^^ = 7 1 TTTT— 77— r ' (^4) 



47fc'(2fc + fc')' 
' + 2{k + k')- 
4 \ mk'^ + {k + k')^ 
1 J mk^ 
4 \mfc2-l- (fc + fc')7 
k 7 

4 [mk'^ +j'^(k + k')] 



1 f my„ . 

^ - . r^.. I'L. , (26) 
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Thus, the contributing terms will be: 

Ci,i, = {7^1(2 + v42)gAA + 7^2(2 + A,}, (27) 

Ci^i2 = ——- {7^1(2 + A2)Qaa+ 7^2(2 + A2)Saa,}, (28) 

z — i(Zi 

Ci,f, = 7^ {7^1(2 + 7^2(2 + A2)}Qaa, (29) 

Ci,v, = —^\ — {2 + Al)nAK + — (2 + Al)ni^A, (30) 



z 






z 


z 


— iCd 




z 


z 


— i(d 




z 


z 


- i(d 




z 


z 


-i(Z) 




z 


z 


-i(Z) 




z 


z 


-i(Z) 



52 = — ^ <! — (2 + ADTZaa + — (2 + ADTZaa } ■ (31) 

Ci,c, = —-H[f,{2 + Al)-f2i2 + Al)] (32) 

z — i 

C^.v, = ■^^-^H[f2i2 + Al)-f,i2 + Al)] (33) 
The other possible terms vanish: 

^Vi V2 — ^Xi Vi — ^X2 V2 — ^' (34) 

since they are integrations of products of odd functions of g's with even functions of g's. 

Notice that the effect of modulation is to scale the temperature value by a factor of 1 + A^/2. This is exactly what 
is obtained by taking the time average of T{1 + Asin{ujt))'^: 



T(l + Asin(wt))2 ^ T{l + A'^/2). 
We notice that the coupling term k is responsible for the non-zero values of Qaa, Qaa, ^aa and 71. 

B. Exact solution for the time-averaged steady state distribution 

Using the results obtained above, it is possible to obtain that Eq. (fT9|) is: 

P-(X1, X2, V,) = rf^^^^.Q,.,+,:P,.,+.Q...+,P... W{Ql,P^,Q2,P2) (35) 

J^oo '^'^ -^"^ ■^'^ 

where: 



OO CXD OO OO 



M=ON=OS=OT=0 

= exp{^H[T,{2 + Al)-T2{2 + Al)] [Q1P2 ~ Q2P1]} x 

X e.J-(9l^AA±lQiQ^^A±QlQ^y,T,(^2 + A^ 



QIQaa + 2O1Q2 Qaa + 0?Qaa 



[7^2(2 + ^2)] 



X exp<; - — 
Pi 

X exp<; - — 



— (2 + ADUaa + — (2 + ADUa a 
m m 

— (2 + ADUaa + — (2 + ADUaa 
m m 



(36) 



The exact final result is given by: 

XiVi +MxiX2XlX2 +Afxi X2'y2 + A/'„l^,2^;l'y2} , (37) 
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where all the coefficients Ma^p^ above depend on the temperatures and mechanical constants of the system. 

Due to the couphngs present in Ea. ([57)) . it does not describe a usual Boltzmann distribution but instead a steady 
state where couplings between position and moments arise. The basic reason for this to occur is that the work done 
by the coupling spring is of the form (work done on particle 1 by the spring) on an interval of time dt 

dWi — —k [xi — X2) vi dt, 

generating, as we shall see, a correlation between X2 and vi due to the coupling above. A similar correlation between 
Xi and V2 also appears. 

It is straightforward to show that these correlation functions are given by: 

(xi«2) = (38) 
{X2vi) = (39) 

where 

^1 — ^ •^XxV2 •^X2X2 -^ViVx J^X\V2 ^X2V\ '^X\X2 •^X2V2 •^V^V-y 

J^XiVi J^X2X2 ^ViV2 •^X\V\ J^X2V\ •I^X2V2 •^V\V2 •^X-iX2 ^X2V\ (^0) 
^2 ■^V2V2 ■^X\X\ ■^X2X2 ■^V\V\ '^'^V\V2 ^X\X\^X2X2 ^^X\X\^X2V2 ^ V\V\ ~^ ^ ^ V\V2 ^ X \X \ ^ X2V\ ^ X2V2 

4 A/ii2f2 ^ X\X\ ^ X2V\ 4 A/^J2^^2 ^X2X2 ^X\V\ ^ ^ X\V2 ^ X2X2 ^V\V\ ~t~ ^^V\V2 ^ X2X2 ^X\V-i •^XiV2 ~t~ 

+ AA ^AA ^ + AA ^AA ^ + 4 AA KT KT KT — 

^ •J^V-iV2 •J^XiX2 ^•^^X2V2 •J^XiVi ^ ^ ^ a;2tJ2 ^ a^i 3:2 ^ a;i 1)2 ^ ^1 ^1 

- 4 AA AA AA "^-IKT \f \f \f + AA ^AA ^- 

2 A/^j^U2 J^X2Vi -^XiX2 '^XiV2 '^•^XiVi •^XiV2 •^X2Vi -^X2V2 ^ AAy2t;2 •^X2V\ •^XiX2 -^XiVi 7 (^-^) 
^3 4 A/a;^a;i ■^X2Vi •^V2V2 ^X2V\ ^X\V2 ^ ^ X\X2 ^X\V\ ^ V2V2 

^ ^X\X\ ^ X2V2 ^ V\V2 ^X2V2 ^X\V-i ^X\V2 •^V\V2 •^X-iX2 •^X-iV2 • (^^) 

Observe that we can swap Pi and P3 by the transformation 1^2. 

In the equilibrium limit Ti ~ T2 — T, all the terms of the form Af^v will vanish. In consequence, the couplings 
represented in Eqs. (j38p and (|39p will also vanish. In equilibrium the flux of heat ceases and velocities decouple from 
positions, as in the following cases shown below. 



C. Interesting limits 

Two interesting limits arise. Firstly, by decoupling the particles 

k = Q ^ Qaa = Qaa = Qaa = 7^AA =n = 0. 
The distribution is given by the product of two independent Boltzmann terms: 

„„, , mk' f fc'x? mv? k'xy mv^) 

P {XUX2,V,,V2) = J-^^ ^''^ 

Secondly, by taking the equilibrium (same temperature) case Ti — T2 — T . The final result corresponds to the 
Boltzmann distribution: 

mWk'(k' + 2k) ( fc'x? k'xn fc(xi — a;^)^ mv? mw? 1 , , 

p^^ix,,X2,v„V2)^ (2.T)^ ^"Pj"^"^" 2T 



V. THERMAL CONDUCTANCE 



We are going to obtain the current of energy (heat) between the two Brownian particles by two methods: the 
exact direct calculation of the work rate between the particles, and a Green-Kubo formalism appropriate for finite 
systems [l8| . 
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The Green-Kubo (GK) formalism '^s'l has apphcations for many problems such as fluid slab flow properties [29| . 
diffusion in granular fluids [10, Ull, |32i], fluctuation-dissipation theory [H, [11], thermal conductance in condensed 
matter systems [1^ [s^, viscosity of trapped Bose gas [s^, triple-point bulk and shear viscosities [H, [1^ or self- 
diffusion [405 for Lennard- Jones fluids, among others. The GK method depends crucially on the convergence of time 
integrations of flux-flux correlation functions. 

The convergence of the GK integral depends on the flux-flux time-correlation functions decaying fast enough, 
otherwise the time integration will diverge such as happens for two dimensional hydrodynamic systems [4l| . This is 
due to the mode-coupling between hydrodynamic modes generating a t~^-dependent tail in the velocities correlation. 
However, for three dimensions the tail goes as t~'^l'^ [HI, [43,[4^ and the Green-Kubo integral converges. On the other 
hand, in one dimension non-diffusive effects can affect the validity of Fourier's law while a Green-Kubo approach 
might still be valid [ll] . 



A. Energy flux 



In order to proceed, we will define the energies and fluxes for our system. The "local" energy density will be given 

by 

£1,2 = ^mw?^2 + ^fc'a;?,2- (45) 

The contact of the particles with the thermal reservoirs, and the presence of the dissipative terms, irnply a flux of 
energy into, and out of, the system at both positions. These instantaneous contact fluxes are given by [l7[: 

Jci = -7«i + (46) 

ic2 = -7U2+V2f72- (47) 

As the coupling spring acts as the interaction channel between the particles, we define, for each particle, the transmitted 
heat flux (Energy/Time) as 

jti = -fc(a;i(<) -X2(t))wi, (48) 
Jt2 = ~k(x^{t)-xx(t))v^. (49) 

The local inter-particle elastic energy is deflned as 

E,i = hz(xr{t)-X2{t)f . (50) 

The total balance of energy requires that the excess energy to be stored in the spring potential. Thus, it is straight- 
forward to see that the above definitions do respect energy balance since 

jti +it2 = -dE^i/dt. 

The effective transfer flux ji2 can now be deflned: 

jl2 = ^(itl - jt2) 

= -Mxi(t)-.2(t))(^^i^^^). (51) 

The deflnition above corresponds to sharing the elastic energy, deflned in Eq. [50l in equal parts between the 
neighboring particles. 



B. Direct calculation of k 



The thermal conductance is: 



^^^(T,AT) = ^(j-i2)at, (52) 
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where Ai = A2 = 0, Ti = T , T2 = T + AT , and ()at is the average at AT > 0. The above expression for k goes 
beyond first order approximation since it contains all the information needed to calculate the heat flux, as shown in 

(ji2)at = (ji2)(T,AT)= / dtKiT,t). (53) 

Jo 

The average heat flux is given by (j'i2}at and can be calculated exactly: 

0'i2)at = -k/ixi - X2) [ ^^^ )/ = (.{X1V2) - {x2Vi)) . (54) 



Using Eqs. ^ and we write: 



k (Pi - P3) 



where the values of (Pi,X'2,2'3) are given in Eqs.lgD]), (gll), and ((i^ . 

After some tedious (but straightforward) algebra, the final result is rather simple: 



2 [to fc2+72(fc + fc')]' 



(^•i2)AT = 2fcHAr^ ^=____^___, (56) 



where k is exact and independent of T and AT. 

It is not unexpected to find the flux proportional to A T, since this result has been obtained for similar models 
before 0, S However, Eg. ([56)1 represents the time-average over the full dynamics of the system. We do not make 
any use of approximate Master-equation- type methods, such as the Fokker-Planck equation to obtain the value 
of K. Our method is equivalent to solving exactly the dynamical equations of motion given the realization of the noise, 
then taking the noise average, and finally time-averaging the final result. In principle, the present approach can be 
generalized for any type of noise, not only white noise. It is interesting to compare Eq. (|56p with the results obtained 
from a Green-Kubo integration. This will be an interesting test on the validity of the choice of the thermal current, 
and also of the approximations used in order to derive the Green-Kubo formalism. 



C. Green-Kubo calculation of k 



The exact expression for k above can be compared with proposals in the literature where Green-Kubo formulations 
for the thermal conductance are given. In the spirit of the prev ious paragraph, the effective flux j plays the role of 
the fluctuating flux ji2 for a Green-Kubo relation proposed for obtaining the thermal conductance: 



lim 



G) 



AT 



1 



at"o at T^ Jo 
where () stands for the equilibrium average (AT = 0). We write 



dt{j{t)jm, 



K = lim — 
O — ^00 iZ 



lim lim 



dt 



1 



(T) 



dr (ji2(i -I- T)ji2{t))AT=o, 



dt e 



dre '''^ (ji2(i + T)ji2(i))AT=o- 



,^0+ 9^0+ (T)2 

Replacing the flux above into Eg. ([55]) . we obtain the Green-Kubo expression for k: 



K = lim lim z 

z^0+ 9^0+ 



dte-'-' — 
2^2 



{Xl{t + t) - X2it + t)) 
[Xi{t) - X2(t)) 



dre X 

Vlit + r) +V2{t + T) 



AT=0 



(57) 



(58) 



(59) 
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After some algebraic manipulation the expression for the thermal conductance becomes: 



~ lim lim lim 



z^o+ e^o+ £^0+ 16 T2 



27r 



27r 



^93 

27r 



(iq-i + e)(i(74 + e) 



z - {{iqx + e) + (i92 + e) + {iq^ + e) + (ig4 + e)] Q- {{iqi + e) + (iga + e) 

X + e) - i2(jgi + e)) (ii(i93 + e) + 52(i'?3 + e))] x 

X [(xi(ig2 + e) - i2(«92 + e)) (5i(ig4 + e) + i2(«g4 + e))]) 

fc^7' r dqi r (^93 + e)("^g3 - 
= hm hm — ■ — / — — / — — r ; -r X 

e^o+ e^o+ 4 27r 2tt 6 - [{iqi + e) + {iq^ + e)] 

1 

+ e) + k][T{-tqi - e) + k][r(tq3 + e) - k][r{-iq3 - e) - k]' 

where the poles are given by: 



(60) 



7 . k' 7^ 
2m V m ' 



S3 



7 . {2k + k') 72 
2m V m 4m2 ' 



and the integration path is shown in FigO 

Equation (pO)) gives exactly the same result of Eq. fSS]) . showing that both approaches are completely consistent. 



/ • i(e-S3) • i(s-s,) 


• i(e-S,') • i(e-S3") \ 






• -q2-q3-q4+3ie 

• -q3-i(e-2e) 


-iz J 



K(q) 



FIG. 2: Integration path over the poles for the Green— Kubo calculation of the conductance. 



D. Discussion 



The coherence shown for the thermal conductance results for a finite systems, calculated either directly, Eg. ([56| . 
or via the Green-Kubo approach, Eq. (j58p . seems to point to the validity of considering the microscopic work as the 
correct fluctuating flux variable to be used for coupled particle systems. In fact, for more realistic models in which 
the number of particles is large, solving the same problem for non-harmonic potentials might be the way to obtain a 
rigorous demonstration of Fourier's Law. 

In our case, despite the somewhat involved aspects of the algebra, the final value for k is quite simple and carries the 
influence of both couplings, k and fc', the friction coefficient 7, and the inertia m. The program we followed in order 
to find K is equivalent to solving the exact equations of motion of the Brownian particles system for each realization 
of the noise functions, and then taking the average over the noise. No approximations of any sort are necessary once 
the basic model is provided. The present treatment can be extended to other distinct kinds of noise, such as colored 
noise (non-Markovian) , or even distinct heath baths acting on the same particles. 

However, the present method can readily be extended to (finite) systems composed of more than two Brownian 
particles, systems that may be large enough to be taken as "macroscopic" . The difficulties to treat such systems are 
operational or numerical, rather than conceptual. 



11 



VI. CONCLUSIONS 



Brownian particles (BP) are an excellent laboratory for studying non-equilibrium physics. They are simple to 

describe but present many of the features of more complex models, such as the possibility of reaching stationary 
states when submitted to thermal contacts at distinct temperature. They are also good approximations for larger 
systems, like polymers, that could be modeled by chains of BP attached to each other by some type of attractive 
potential. 

Another interesting characteristic of such systems is that they are simple enough so that we can extract exact 
solutions for their long-time behavior. This allows us to obtain results that are hard to come by using other methods. 
It is already known that we can obtain exactly the equilibrium probability distribution for Brownian particles subjected 
to Markovian or non-Markovian noise, or a combination of both. This type of external forcing allows us to keep a 
system formed by a single particle constantly on an out of equilibrium steady-state. 

Furthermore, techniques based on time-averaging are very interesting since they are ensemble independent, driven 
only by the dynamical relations governing the interaction Brownian particle-heath bath. In fact, this corresponds to 
following a system during the realization of an experiment. 

In the present work, we have studied the thermal conductance for a system of coupled particles, by taking advantage 
of the mechanically simple characteristics of Brownian particles and of time-averaging. Our system consists of two 
particles coupled by a spring potential, with the heat flux flow 712 being due to the mechanical work done through the 
spring coupling the two particles. The particles may be kept in contact with distinct thermal baths and the flow of 
energy may be obtained by means of the thermal conductance coefficient k, appropriate for small systems. The latter 
is calculated both from first principles, and by means of a Green-Kubo formulation, and the obtained exact results 
are identical. The final form of k depends only on the variables of the system, such as the mass of the particles, the 
spring couplings and the friction coefficients. The thermal conductance can be thought as a first step for obtaining 
the equivalent form for, more sophisticated, macroscopic systems such as long polymers. 

We believe that the coherence between the exact direct calculation and the exact Green-Kubo formulation for k 
shows the correctness of the basic definitions, in special that of the heat flux, used in the problem. 
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